Double gamma-ray lines from unassociated Fermi-LAT sources revisited 
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We search for the presence of double gamma-ray line from unassociated Fermi-LAT sources in- 
cluding detailed Monte Carlo simulations to study its global statistical significance. Applying the Su 
& Finkbeiner selection criteria for high-energy photons we obtain a similar excess over the power- 
law background from 12 unassociated sources. However, the Fermi-LAT energy resolution and the 
present low statistics does not allow to distinguish a double peak from a single one with any mean- 
ingful statistical significance. We study the statistical significance of the fit to data with Monte 
Carlo simulations and show that the fit agrees almost perfectly with the expectations from random 
scan over the sky. We conclude that the claimed high-energy gamma-ray excess over the power-law 
background from unassociated sources is nothing but an artifact of the applied selection criteria and 
no preference to any excess can be claimed with the present statistics. 



Introduction. The evidence for 130 GeV gamma-ray 
excess from Galactic centre consistent with dark matter 
(DM) annihilations in the central cusp is by now well es- 
tablished fact [TH3]. Therefore searches for the signal of 
DM annihilations from other DM dominated objects is of 
utmost importance. An evidence for the same 130 GeV 
peak has been found from nearby galaxy clusters support- 
ing the DM annihilations as an origin of the peak [4] . The 
two signals come from known DM dominated sources. A 
search for possible unidentified sources of the 130 GeV 
gamma-ray peak was performed in [3J where several of 
them have been identified. It is attempting to associate 
this signal with DM subhaloes of our Galaxy. However, 
due to large "look elsewhere effect" related to uniden- 
tified sources those can just be statistical upward fluc- 
tuations of the background [BJ. More data is needed to 
resolve this issue. 

Recently a similar search for possible gamma-ray ex- 
cess from unassociated Fermi-LAT point-like sources was 
performed in [BJ. The motivation for that search was a 
claim that possibly those sources can be identified with 
DM subhaloes, exactly the same motivation as in [2J. 
The smart trick the authors use is to search for a double 
peak [7] motivated by most of particle physics models 
of DM annihilations to photons. If such a double peak 
is observed, this would favour particle physics origin of 
the excess over astrophysics. Although the present low 
Fermi-LAT statistics together with its limited energy res- 
olution does not allow to distinguish a double peak struc- 
ture from a single broad peak with a meaningful statis- 
tical significance [S], the double peak gives marginally 
better fit to Galactic centre data [3]. The authors of [BJ 
claim to have observed a double peak with 3.3<r local 
statistical significance over the power-law background. 

The analyses and results of [BJ can be criticized based 
on two arguments. First, the Fermi unassociated sources 
are identified by gamma-rays with much lower energy 



than the 130 GeV peak. Consequently the first ques- 
tion to ask is why those sources play any role at high- 
energies and what is the possible connection between the 
low and high-energy gamma-rays. Those issues have been 
addressed in [91 [TO] . Second, the authors of [BJ make no 
attempt to estimate the global significance of their claim. 
While stacking together data from several regions, they 
should have estimated the expected effect of their selec- 
tion procedure by choosing arbitrarily the same number 
of regions in the sky, and scan over the full area as is 
done in [3]. This is the issue addressed in this note. 

The aim of this work is first to repeat the analyses of 
6J using their criteria for selecting high-energy photons, 
and then to apply the Monte Carlo simulations similar 
to the one performed in [3] to study their claim. We 
do find a similar excess of high-energy gamma-rays over 
the power-law background as in [BJ with a similar local 
significance confirming that the idea of looking for the 
double peak over the power-law background is indeed 
very useful one. However, we show that with Fermi-LAT 
energy resolution and with present very limited statis- 
tics one cannot distinguish between the double and sin- 
gle peaks with any meaningful significance. We do not 
observe any systematic effects in the signal that could be 
coursed by the detector since those must have occurred 
everywhere in the data, including background. Thus the 
claims of detector effects in [HI [TO] have no justification at 
present. Performing the Monte Carlo study to determine 
the global significance of the excess we find that the se- 
lection criterion of [BJ itself implies a peak, and the peak 
from the unassociated Fermi sources agrees almost per- 
fectly with the Monte Carlo expectation. We conclude 
that the claimed excess is nothing but an artifact due to 
the applied event selection and is not associated with the 
Fermi-LAT point-like sources. Thus the global statistical 
significance of the excess is completely negligible and the 
results of [BJ are systematically biased. 
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Data analyses and Monte Carlo results. The 

authors of [BJ work with 319 unassociated Fermi-LAT 
sources and select a subset of them (16 sources) which 
have at least one photon in the energy range 100- 
140 GeV. We find 12 objects like that using the data 
selection constraints as described in [5] and suggested by 
Fermi-LAT team, which is 4 less than found in [BJ (but 
in agreement with [9]). We note that the total number of 
photons in the energy range 20-300 GeV from those 12 
sources is 37 that is really small number for any meaning- 
ful statistics. For comparison, in our work [2J we did not 
consider any individual sources with less than 80 high- 
energy photons at all. Nevertheless, stacking together 
photons from the 12 sources we observe an excess over 
power-law background depicted with continuos red line 
in Fig. [T] This is very similar to the result obtained in 
[BJ. Closer look at these 12 sources reveal that there are 
only 12 photons in the energy range 100-140 GeV. In- 
terestingly, 9 of them are either close to the 111 GeV 
(108.9-116.6 GeV) or 129 GeV (124.7-133.4 GeV) en- 
ergy intervals. We estimate that the probability of this 
happening is rather low (1.5%). This result is the reason 
behind the claims of [BJ observing the double line with 
high local statistical significance. 

However, notice that the fit of data in Fig. [I] (the red 
continuous line) does not posses a double peak. The rea- 
son is that physically meaningful kernel size for this anal- 
yses is determined by the Fermi-LAT energy resolution. 
We remind for the reader that our analyses is based on 
kernel smoothing method that is independent of binning 
and refer the reader to Ref. [2] for technical details. If 
we choose narrower kernel we can obtain the double peak 
structure as depicted in Fig. [T] with red dashed line. How- 
ever, such a choice is not physically well motivated at 
the moment because of limited Fermi-LAT energy reso- 
lution. Playing with the existing data (choosing small 
kernel sizes in our case or choosing binning in the case 
of [BJ) one can, indeed, see some evidence for the double 
peak. However, its statistical significance over a singe 
peak is just marginally better and the present Fermi- 
LAT energy resolution together with low statistics does 
not allow to distinguish the double peak from a single 
one. This criticism is completely general and applicable 
to all searches for the double peak in Fermi-LAT data. 
Only the future experiments with high energy resolution 
and high statistics can really tell the two peaks apart. 

Coming to the main motivation of this work, a selection 
of 12 unassociated sources (that are identified by low- 
energy gamma rays) with one energetic photon cannot 
be physically justified. To the contrary we expect that 
such a procedure is just choosing the expected signal out 
of 319 samples by hand. Repeating the same with any 
12 energetic photons in the sky should give the similar 
peak. Thus we expect such a "signal" to be an artifact of 
the selection of photons and not to be a signal associated 
with the particular objects. The authors of [BJ do not 



— Monte Carlo average (95% CL) 

Fermi 12 unassoc. sources (smooth) 

- Fermi 12 unassoc. sources (detail) 

Background (2.6 power law) 





















30 100 130 200 

E [GeV] 



FIG. 1: Fit to Fermi-LAT gamma ray spectrum from 12 unas- 
sociated sources selected as described in the text (solid red 
curve) together with the average Monte Carlo expectation 
from 12 arbitrary regions in the sky selected by the same cri- 
teria (solid green curve). The 2a error- band obtained with 
Monte Carlo technique is presented by the grey band. Fit to 
data with unphysical energy resolution (narrow kernel) is also 
presented with dashed red line. The excess claimed in Ref. 6 
is clearly an artifact of the selection criteria of photons and 
can be obtained for arbitrary selection in the sky. 



make an attempt to estimate statistical significance of 
the excess using Monte Carlo methods by scanning over 
the sky as we did in [4] . Instead they limit themselves to 
the 16 point-like sources. This is not a valid statistical 
sample and cannot be used for estimating the statistical 
significance of the excess. 

To show that our criticism is well justified we use the 
selection criteria of [BJ and repeat their analyses for 12 
regions in the sky using Monte Carlo method described 
in [3]. The results are presented in Fig. [I] The green 
line is an average of Monte Carlo result, this is what one 
typically expects to obtain if one repeats the procedure 
of [BJ for any 12 arbitrary regions in the sky. This agrees 
surprisingly well with the fit to data of 12 unassociated 
sources. In fact, by comparison of the red and green 
lines one can conclude that the data from 12 unassoci- 
ated sources follows the Monte Carlo expectations almost 
perfectly. The 2a error band of the Monte Carlo result 
is presented in Fig. [T] with the grey band. The data is 
well within the band. We conclude that the observed ex- 
cess has nothing to do with the point like sources, this 
is just an artifact of the procedure of choosing the pho- 
tons. Therefore the statistical significance of the excess 
claimed in [BJ is clearly and systematically overestimated. 
In fact, our results show that no excess can be claimed 
at all. 

For completeness we also studied low-energy gamma- 
ray spectra of the 12 unassociated sources and find their 
spectral shapes to be all very different from each other. 
If the low energy photon spectrum in those regions is 
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dominated by the continuum spectrum of the Z'y, h'y 
final states of DM annihilations, the low-energy spectra 
must also be identical. We conclude, in agreement with 
[5] [TU] , that this is not the case. 

We finally comment that the 16 candidate sources 
found in [5J are mostly located close to the galactic disc, 
exactly as the sources in [2] . Whether those are DM sub- 
haloes or not requires independent checks. The present 
poor statistics does not to allow for definite conclusion. 

Discussion and conclusions. The claim that high- 
energy photons from Fermi unassociated sources orig- 
inate from DM annihilations in DM subhaloes of our 
Galaxy has been criticized in [9] [10] based on the shape 
and magnitude of the low energy photon signal from the 
selected 12 unidentified sources. Both those papers show 
that at most 2 sources out of 12 can be potentially asso- 
ciated with DM subhaloes, and most probably the low- 
energy signal comes from AGNs [10] . While this criticism 
might be relevant, it is based on the model dependent as- 
sumption that the high-energy and low-energy photons 
are physically associated. For example, it is logically pos- 
sible that the Fermi unassociated sources are, indeed, 
misidentified astronomical sources which are within DM 
subhaloes. 

To explain the findings of [6] the authors of [9] [10] have 
also argued that the double peak in unassociated sources 
is potentially due to some kind of unknown artifact of 
Fermi-LAT detector without showing what this artifact 
could be or how it would give the double peak. Moreover, 
from Monte Carlo simulations (from arbitrary regions in 
the sky) we do not observe this double peak, showing 
that the double peak is not visible everywhere as would 
be expected if this is unknown artifact of Fermi-LAT de- 
tector. 

Since the low energy data cannot rule out the possi- 
bility that high-energy signal comes from DM subhaloes, 
we have studied the claims related to the high-energy 
photon signal in detail. We find that choosing the 12 
Fermi unassociated sources using the selection criteria for 
high-energy gamma-rays as in [6], stacked data, indeed, 
shows a peak over the power-law background. For the 12 
sources there are only 12 photons in the energy range 100- 
140 GeV, and 9 of them are close to either the 111 GeV 
or the 129 GeV energies. Probability for this to happen 
is less than 2% (according to Monte Carlo simulations 
considering only the high-energy range 100-140 GeV), 
and this is the reason why [5] claim the double peak in 
data. However, we show that the energy resolution of 
Fermi LAT, together with poor statistics, does not allow 
to to distinguish the double peak from one peak with 
any meaningful significance, as demonstrated in Fig. [T] 
At present this result applies to any search for double 
peaks in Fermi LAT data. 

The most important result of this work is that we have 
demonstrated with Monte Carlo simulations that the ex- 
cess from the 12 unassociated Fermi-LAT sources is in 



a perfect agreement with expectations if one scans ran- 
domly over sky and chooses 12 energetic photons accord- 
ing to the criteria of [6] . Thus the excess over power-law 
background is an artifact of selecting high-energy pho- 
tons out of the signal sample and has nothing to do 
with the unassociated sources. Therefore the choice of 
16 point-like sources as an statistical ensemble, as done 
in [6] , greatly overestimates the statistical significance of 
the excess. Since the fit to data agrees with the Monte 
Carlo expectation well within the statistical error band 
presented in Fig. [I] no claim of any excess can be made 
with present statistics. 

We also searched for possible detector effects that could 
be the source of the double line, as speculated in [SJ [TU] , 
and do not find any. If there is a systematic detector 
effect, it must show up not only in the 12 chosen regions 
but also in the background. There is no such an effect 
in the statistically much larger background. We conclude 
that most probably the observed double structure of the 9 
high-energy photons is just a statistical fluctuation. More 
data from new experiments with high energy resolution 
should solve the open issue of one peak versus double 
peak in the observed 130 GeV excess from the Galactic 
centre and from nearby galaxy clusters. 

We stress that the results and arguments presented in 
this note do not affect any way the evidence for 130 GeV 
excess from Galactic centre and from nearby galaxy clus- 
ters. 

Acknowledgement. We thank Doug Finkbeiner 
for useful communication. This work was supported 
by the ESF grants 8090, 8499, 8943, MTT8, MTT59, 
MTT60, MJD52, MJD272, by the recurrent financing 
projects SF0690030s09, SF0060067s08 and by the Euro- 
pean Union through the European Regional Development 
Fund. 



* andi.hektor@cern.ch 

T martti.raidal@cern. ch 

* elmo@aai.ee] 

[1] C. Weniger, |arXiv:1204.2797| [hep-ph] . 
[2] E. Tempel, A. Hektor and M. Raidal, |arXiv:1205.l"04"5] 
[hep-ph]. 

[3] M. Su and D. P. Finkbeiner, |arXiv:1206."T616| [astro- 
ph.HE]. 

[4] A. Hektor, M. Raidal and E. Tempel, |arXiv:1207.4466| 
[astro-ph.HE]. 

[5] A. Boyarsky, D. Malyshev and O. Ruchayskiy, 
|arXiv:1205.4700] [astro-ph.HE], 

[6] M. Su and D. P. Finkbeiner, |arXiv:1207.7060| [astro- 
ph.HE]. 

[7] J. M. Cline, |arXiv:1205.2688| [hep-ph], 

[8] A. Rajaraman T. M. P. Tait and D. Whiteson, 

|arXiv:1205.4723] [hep-ph]. 
[9] D. Hooper and T. Linden, |arXiv:1208.0828| [hep-ph], 
[10] N. Mirabal, [arXiv:1208.1693| [astro-ph.HETT 



